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ABSTRACT 

We use high-quality optical rotation curves of 9 low-luminosity disk galaxies to obtain 
the velocity profiles of the surrounding dark matter halos. We find that they increase 
linearly with radius at least out to the edge of the stellar disk, implying that, over the 
entire stellar region, the density of the dark halo is about constant. 

The properties of the mass structure of these halos are similar to those found for 
a number of dwarf and low surface brightness galaxies, but provide a more substantial 
evidence of the discrepancy between the halo mass distribution predicted in the Cold 
Dark Matter scenario and those actually detected around galaxies. We find that the 
density law proposed by Burkert (1995) reproduces the halo rotation curves, with halo 
central densities {po ^ 1-4 xlO~^^ g cm~'^) and core radii (ro ^ 5-15 kpc) scaling as 

-2/3 
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1 INTRODUCTION 

Rotation curves (RCs) of disk galaxies are the best probe 
for dark matter (DM) on galactic scale. Notwithstanding 
the impressive amount of knowledge gathered in the past 20 
years, some crucial aspects of the mass distribution remain 
unclear. In fact, the actual density profile of dark halos is 
still matter of debate; we do not even know whether it is 
universal or related to some galaxy property, such as the 
total mass. This is partly because such issues are intrin- 
sically crucial and partly because it is often believed that 
a RC leads to a quite ambiguous information on the dark 
halo density distribution (e.g. van Albada et al. 1985). How- 
ever, this argument is true only for rotation curves of low 
spatial resolution, i.e. with no more than 2 measures per ex- 
ponential disk length-scale Rd, as is the case of the great 
majority of HI RCs. This is because the galaxy structure 
parameters are very sensitive to both the amplitude and 
the shape of the rotation curve in the region < r < Rn, 
(e.g. Blais-Ouellette et al., 1999) which corresponds to the 
region of the RC steepest rise. No reliable mass model can 
be derived if such a region is poorly sampled and/or radio 
beam-biased. However, in case of high-quality optical RCs 
with tens of independent measurements in the critical re- 
gion, the kinematics can probe the halo mass distribution 
and resolve their structure. 

Since the dark component can be better traced when 
the disk contributes to the dynamics in a modest way, it is 
convenient to investigate DM-dominated objects, like dwarf 
and low surface brightness (LSB) galaxies. It is well known 
that for the latter there are strong claims of dark matter 
distributions with cores of constant density, in disagreement 
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116-G12 
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-20.0 


531-G22 
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-20.4 


563-G14 


2.0 


148.9 


-20.5 


79-G14 


3.9 


167.1 


-21.4 


M-3-1042 


1.5 


148.0 


-20.1 


N7339 


1.5 


172.7 


-20.6 


N755 


1.5 


102.4 


-20.1 



Table 1. Observational properties of the sample galaxies: (1) 
Galaxy name. (2) Length— scale of the exponential thin disk (from 
Persic & Salucci 1995). The optical radius Ropt is the radius en- 
compassing 83% of the light: Ropt = 3.2Riy. (3) Circular velocity 
at Ropt- (4) Total /-band absolute magnitude (Mathewson, Ford 
& Buchhorn, 1992). 



with the steeply cusped density distributions of the Cold 
Dark Matter Scenario. (Flores & Primack, 1994; Moore, 
1994; Burkert, 1995; Burkert & Silk, 1997; Kravtsov et al., 
1998; McGaugh & de Blok, 1998; Stil, 1999). However, these 
findings are 1 ) under the caveat that the low spatial resolu- 
tion of the analysed RCs does not bias the mass modeling 
and 2) uncertain, due to the limited amount of available 
kinematical data (see van den Bosch et al., 1999). 

In this paper we will investigate the above-discussed 
issue by analysing a number of high-quality optical rota- 
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Figure 1. Rotation curves of the sample galaxies. 



tion curves of low luminosity late-type spirals taken from 
Persic and Salucci (1995, PS95), with I-band absolute mag- 
nitudes —21.4 < Mi < —20.0 that, in terms of rotational 
velocities, translates into 100 < Vopt < 170 km s~^. Objects 
in this luminosity /velocity range are DM dominated (e.g. 
Persic, Salucci & Stel, 1996), but their RCs, measured at 
the PS95 angular resolution of 2", have a spatial resolution 
of w ~ 100(/)/10 Mpc) pc and ridata ~ Ropt/w indepen- 
dent measurements. For nearby galaxies: w « Rd and 
ndata > 25. Moreover, we select rotation curves of bulge- 
less systems, so that the stellar disk is the only baryonic 
component lor r ^ Rd- 

Since most of the properties of cosmological halos are 
claimed universal, it is worth to concentrate on a small and 
particular sample of RCs, that however can provide crucial 
information on the dark halo density distribution. The sys- 
tematics and the cosmic variance of the DM halos are inves- 
tigated elsewhere (Salucci and Burkert, 1999; Salucci, 1999). 
Finally, let us stress that to estabilish the actual theoreti- 
cal properties of CDM halos or to investigate non standard 
CDM scenarios possibly in agreement with observations is 
beyond the scope of this work. 

In §2 we describe our sample of RCs. We present our 
mass modeling technique and its results in §3. In §4 we dis- 
cuss the inferred halo profiles and their properties. In §5 
wo perform a disk-halo rotation curve modeling by adopt- 
ing a CDM density profile for the dark halo. We summarize 
our results in §6. Throughout this paper we adopt a Hubble 
constant of Hq = 75 km s~^ Mpc~^. 



2 ROTATION CURVES 

The rotation curves of the Persic & Salucci (PS95) 'excel- 
lent' subsample of 80 galEixies are all suitable for an accurate 
mass modeling. In fact, these RCs properly trace the gravi- 
tational potential in that: 1) data extend at least to the op- 
tical radius, 2) they are smooth and symmetric, 3) they have 
small rms, 4 ) then' have high spatial resolution and a homo- 
geneous radial data coverage, i.e. about 30 — 100 data points 
homogeneously distributed with radius and between the two 
arms. From this subsample we extract 9 rotation curves of 
low luminosity galaxies (5 x Lq < Li < 2 x 10^" L0; 
100 < Vopt < 170 km s~^), with their /-band surface lumi- 
nosity being an (almost) perfect radial exponential. These 
two last criteria, not indispensable to perform the mass de- 
composition, are however required to minimize the uncer- 
tainties of the resulting dark halo density distribution. The 
selected RCs are shown in Figure 1 (for all details we refer to 
PS95). They are still growing at Ropt, in that mostly tracing 
the dark halo component. Each RC has 7—15 velocity points 
inside Ropt, each one being the average of 2 — 6 independent 
data. The RC spatial resolution is better than 1/20 Ropt, the 
velocity rms is about 3% and the RCs logarithmic derivative 
is generally known within about 0.05. 

3 MASS MODELING 

We model the mass distribution as the sum of two compo- 
nents: a stellar disk and a spherical dark halo. By assuming 
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tial thin disk. Light traces the mass via an assumed radially 

constant mass-to light ratio. 

In the r.h.s of cq.(l) we neglect the gas contribu- 
tion Vgas{r) since in normal spirals it is usually modest 
within the optical region (Rhee (1996), Fig. 4.13): figaa = 
(^gas/^^)i?opt ~ 0.1. Furthermore, high resolution HI obser- 
vations show that in galaxies with RCs similar to those of the 
present sample ( M33: CorbeUi & Salucci, 1999; NGC300: 
Puche, Carignan & Bosma, 1990; N5585: Cote & Carignan, 
1991; N3949 and N3917: Verheijen, 1997) Vga^r) is well rep- 
resented by Vgas{r) ~ for r < Rd and: 



Vgas{r) ~ (20 ± 5)(r - Rd)/2Rd 



Rd <r< 3Rd{2) 



Table 2. CDRmass models: (1) Galaxy name; (2)-(3) Parameters 
of the best-fit CDR models with their la uncertainties (§3); (4) 
Disk mass-to-light ratio in the /-band in solar units. 



centrifugal equilibrium under the action of the gravitational 
potential, the observed circular velocity can be split into 
these two components: 



V^ir)^VE{r) + V^{r) 



(1) 



This is true as long as one assumes the galactic disk be 
in dynamic equilibrium. By selection, the objects are bulge- 
less and the stellar component is distributed like an exponen- 



Thus, in the optical region: i) Vgas{r) << V^{r) and 
li) d{V'^{r) — Vga,,{r))/dr ^ 0. This last condition implies 

that by including Vgas in the r.h.s. of cq.(l) the halo veloc- 
ity profiles would result steeper and then the core radius in 
the halo density larger. Incidentally, this is not the case for 
dwarfs and LSBs: most of their kinematics is affected by the 
HI disk gravitational pull in such a way that neglecting it 
could bias the determination of the DM density. 

The circular velocity profile of the disk is (Freeman, 
1970): 



p2 
^opt 



{IqKo - /lJ^l)l.6r-/ilopt 
{loKo - IlKl)l.6 



(3) 
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Figure 3. Relation between the halo parameters, (o is in units 

of Rovt) 



where /„ an are the modified Bessel functions, Vopt 
is the measured circular velocity at Ropt and /? = ( ^ J 

^ ^ Ropt 

The parameter (3 represents the disk contribution to the 
total circular velocity at Ropt and can vary from to 1. On 
grounds of simplicity, we choose /3 as the disk free parameter 
rather than the disk mass to light ratio. 

For the DM halo we assume a spherical distribution, 
whose contribution to the circular velocity Vfr(r) is given 
by (Persic, Salucci & Stel 1996, Salucci 1997): 



(a;2 + a2) 



(4) 



whore x = r / Ropt and a is the core radius measured in 
units of Ropt- From eq.(l): 7 = (1 — 13). Since we normalize 
(at Ropt) the velocity model (V^ + V^)^^^ to the observed 
rotation speed Vopt, P enters explicitly in the halo veloc- 
ity model and this reduces the free parameters of the mass 
model to two. 

It is important to remark that, out to Ropt, the proposed 
Constant Density Region (CDR) model of eq.(4) is neutral 
with respect to the competing mass models. Indeed, by vary- 
ing P and a, it approximately reproduces the maximum- 
disk, the solid-body, the no-halo, the all-halo, the CDM 
and the coreless-halo models. For instance, CDM halos with 
concentration parameter c = 5 and Vs = Ropt (see §5) are 
well fit by cq.(4) with a ~ 0.33. 

For each galaxy, wc determine the values of the param- 
eters P and a by means of a x'^~iiiinimization fit to the ob- 
served rotation curves: 



Vmodei{r;P,a) 



V^(r;/3) + V|(r;/3,a) 



(5) 



A central role in discriminating among the different 
mass decompositions is played by the derivative of the ve- 
locity field dV/dr. It heis been shown (e.g. Persic & Salucci 
1990b, Persic & Salucci 1992) that by taking into account 
the logarithmic gradient of the circular velocity field defined 
as: 



V(r) 



dlogVjr) 
rflogr 



(6) 
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Figure 4. N755 rotation curve compared with over-maximum- 
disk models (/3 = 0.2 and a = 2.0 {top) or a = 1.0 (bottom)). Solid 
line represents the total circular velocity, short- and long— dashed 
lines represent the disk and halo contributions. The best-fit model 
for this galaxy has B = O OQ^q JJg and a = l.O^g'j, corresponding 
to Mjj/Lz = 0.2 Mq/Lz 0. We realize that models with signif- 
icantly higher disk mass-to-light ratio are inconsistent with the 
inner rotation curve, regardless of the halo profile. 



one can significantly increase the amount of informar 
tion available from kinematics and stored in the shape of 

the rotation curve. So wc consider x^-s calculated on both 
velocities and logarithmic gradients: 



2 

Xv 



E 

i=l 



Vi-Vm 



idel 



(.n;l3, a) 



5Vi 



V(r'j) - Vmodei{ri\ (3, a) 
5Vi 



(7) 



(8) 



where Vmodei {ri,(3, a) is computed from eq. (3) — (4) and 
eq.(6). As the linear combination of x^-s still follows the 

X^-statistics (Bovingtou & Robinson 1992), we derived the 
parameters of the mass models by minimizing a total Xtot; 
defined as: 



2 _ 2 , 2 
Xtot = Xv + Xv 



(9) 



Xtot 



The above is the a = 1 case of the general relation 
= Xv + • Xv w6 adopt in that the circular velocity 
at radius Vi and the corresponding log- gradient V(ri) are 
statistically independent. Notice that the Xtot best-fit solu- 
tions are not statistically difi^erent from those obtained with 
the usual x^ procedure, however, the ellipse uncertainty is 
now remarkably reduced. 

The parameters of the best-fit models are listed in Table 
2, along with their Icr uncertainties. The mass models are 
well specified for each object: the allowed values for (3 and 
a span a small and continuous region of the (a, (3) space. 
We get a "lowest" and a "highest" halo velocity curve by 
subtracting from V(r) the maximum and the minimum disk 
contributions Voir) obtained by substituting in eq.(3) the 
parameter (3 with f3best + Sf3 and Pbest — 5(3, respectively. 
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Figure 5. (left): The halo velocity profiles of the sample galaxies. V^f (r) rises almost linearly with radius: the DM halo density remains 
approximately constant (right): halo central density po vs. core radius ro- The dashed line is the Burkert relationship for dwarf galaxies. 



The derived mass models are shown in Figure 2, along- 
side with the separate disk and halo contributions. It is then 
obvious that the halo curve is steadily increasing, almost lin- 
early, out to the last data point. The disk-contribution f3 and 
the halo core radius a span a range from 0.1 to 0.5 and from 
0.8 to 2.5, respectively. In each object the uniqueness of the 
resulting halo velocity model can be realized by the fact that 
the maximum-disk and minimum-disk models almost coin- 
cide. In Figure 3 we show the correlation between the halo 
parameters: halos which axe more dynamically important at 
Ropt (i e. with higher 1 — /3) have smaller core radii. Part 
of the scatter of the relation may arise because the sample 
is limited in statistics and in luminosity range (see Salucci 
and Burkert, 1999). Remarkably, we find that the size of the 
halo density core is always greater than the disk characteris- 
tic scale length Rd and it can extend beyond the disk edge 
(and the region investigated). 

As regards to the HI disk contribution, we checked the 
consequence of neglecting it. Its contribution to the circular 
velocitiy (see §3, Corbelli and Salucci, 1999) is computed and 
then actually subtracted from V{r). As result, the values of 
the parameter (3 do not change whereas we obtain slightly 
larger values of the core radius a. As a typical example: for 
N755, considering the gas distribution we get: (3 = O.OS^QlQg 
and a = l.llo^2 to be compared with the value in Table 2. 

Although the present sample is small for a thor- 
ough investigation of the stellar mass-to-light ratios, it 
is worth noticing that the disk mass-to-light ratios found 
(< Md/Li >~ Q.9>Mq / Liq) are typical of late-type spirals 
with young dominant stellar populations and ongoing star 
formation (e.g. de Jong, 1996), which adds to variations in 
stellar populations due to difiierences in age, metallicity, star 
formation history and to uncertainties in the estimate of dis- 
tances and/or internal extinctions In detail, wo ensured that 
for the three galaxies with the lowest disk mass-to-light ra- 
tios, a significantly larger value of /3 was inconsistent with 
the (inner) rotation curve (see Figure 4). 



Finally, it is worth stressing that the present analysis 
computes the mass-to-light ratio from the model param- 
eters as a secondary quantity; the uncertainties in Col. (4) 
of Table 2 do not indicate the goodness of the halo mass 
models, but only specify how well we know this quantity. 



4 DARK HALOS PROPERTIES 

In Figure 5 {left) we show the halo velocity profiles for the 
nine galaxies. The halo circular velocities are normalized to 
their values at Ropt and expressed as a function of the nor- 
malized radius r/Ropt- These normalizations allow a mean- 
ingful comparison between halos of difi'erent masses: the ra- 
dius scaling removes the intrinsic dependence of size on mass 
(more massive halos are bigger) , whereas the velocity scaling 
takes into account that more luminous galaxies have higher 
circular velocities. It is then evident that the halo circular 
velocity, in every galax;y, rises almost linearly with radius, 
at least out to the disk edge: 

VH{r) DC r 0.05i?opt <r< Ropt (10) 

The halo density profile has a well defined core ra- 
dius within which the density is approximately constant. 
This is inconsistent with the singular halo density distri- 
bution emerging in the Cold Dark Matter (CDM) scenario 
of halo formation (Navarro, Frenk & White (NFW) 1995, 
1996, 1997; Cole & Lacey, 1997; Tormen et al., 1997; Tissera 
& Dominguez-Tenreiro, 1998; Nusser & Sheth, 1999). More 
precisely, since the CDM halos are, at small radii, likely more 
cuspy than the NFW profile: pcDM oc r'^'^ (Fukushige & 
Makino, 1997; Moore et al, 1998; Jing, 1999; Jing & Suto, 
1999; Ghigna et al., 1999), the steepest CDM halo veloc- 
ity profile yff(r) cx r^^'^ results too shallow with respect to 
obervations and therefore inconsistent with eq.(lO). Burk- 
ert (1995) has proposed for the halo density distribution the 
following phenomenological profile: 
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Figure 6. NFW best-fits solid lines of the rotation curves (filled circles) compared with the CDR fits (dashed lines). The values are 
also indicated. 



ps(r) = 



(11) 



then, in order to not violate the constraint given by eq.(lO), 
all of them should have the same solid body velocity profile. 



where po (the central density) and ro (the scale radius) are 
free parameters. This density law has a core radius of size 
ro, and, at largo radii, converges to the NFW profile. The 
two parameters po and ro are found correlated: po ~ r^^^^, 
so that the halo profiles reduce to a one-parameter family 
of curves (Burkert, 1995). 

We find that, by adjusting po and ro, we can reproduce, 
over the available radial range, the halo velocity profiles of 
eq.(4) and Table 1. We show in Figure 5 (right) the pairs 
(po, ro) corresponding to the 9 objects. The typical uncer- 
tainties on Po and ro are estimated at a level of about 0.15 
dex and 0.07 dex, respectively: we then confirm the struc- 
ture relation Burkert (1995) found for 5 dwarf galaxies (see 
Figure 5 right). 

Let us stress that the CDR halo properties of eq.(lO) 
raise important issues by themselves and make quite irrele- 
vant, in comparing theory and observations, arguments per 
se important such as the CDM halos cosmic variance, the 
actual value of the concentration parameter or the effects 
of baryonic infalls or outflows (e.g. Martin, 1999; Gelato & 
Sommer-Larsen, 1999). 

We also want to remark that our finding implies that 
disk galaxies are embedded (inside Ropt) in a single dark 
halo. Indeed, if more dark halos were relevantly present, 



5 THE INADEQUACY OF CDM MASS 
MODELS 

Although the mass models of eq.(4) converge to a distribu- 
tion with an inner core rather than with a central spike, 
i.e. to a ~ 1 rather than to a ~ 1/3, it is worth, given 
the importance of such result, to also check in a direct way 
the (in)compatibility of the CDM models with the observed 
kinematics. We assume the two-parameters functional form 
for the halo density by Navarro, Frenk & White (NFW, 1995, 
1996, 1997): 



PNFw(r) = 



(r/r,)(l + r/r,)2 



(12) 



where Vs is a characteristic inner radius and ps the cor- 
responding density. In order to translate the density pro- 
file into a circular velocity curve for the halo, we make use 
of virial parameters: the halo virial radius Rvir, defined as 
the radius within which the mean density is Avir times the 
mean universal density pm at that redshift, and the associ- 
ated virial mass Mvir and velocity Vvir = GMvir/RviT. By 
defining the concentration parameter as Cvir = Rvir/rs the 
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Mvir < 2 • 1012 M0 
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Table 3. NFW mass models: (1) Galaxy name; (2)-(3)-(4)-(5) Parameters of the NFW best-fit models with their 1<t uncertainties. 
An upper limit of 2 • W^'^Mq is imposed on M^ir- Disk mass-to-light ratios in the /-band in Col. (5) are in solar units; (6)-(7)-(8)-(9) 
Parameters of the best-fit models with no constraint on the virial mass. 



halo circular velocity VcY>m{r) takes the form (Bullock et al. 
1999): 



A{x) 



'^(Cvir) 



(13) 



As the 



where X = r/va and A{x) = ln(l -|- x 
relation between Vvir and Rvir is fully specified by the back- 
ground cosmology, we assume the currently popular ACDM 
cosmological model, with Clm ~ 0.3, f^A = 0.7 and h = 0.75, 
in order to reduce from three to two (cvir and r^) the inde- 
pendent parameters characterizing the model. According to 
this model, Avir — 340 at « ~ 0. The choice is conservative: 
a high density Qm = 1 model, with a concentration parame- 
ter Cvir > 12, is definitely unable to account for the observed 
galaxy kinematics (Moore 1994). 

Though N-body simulations and semi-analytic investi- 
gations indicate that the two parameters Cvir and rs seem to 
correlate, we leave them independent to increeise the chance 
of a good fit. Since the objects in our sample are of low lumi- 
nosity, i.e. L,/12 ^ Li < Z/,/3, we conservately set for the 

halo mass Myir an upper limit of Mup — 2 x 10^^ M©, com- 
parable to the total mass of the Galaxy (e.g. Wilkinson & 
Evans, 1999) and to the mass of bright spirals in pairs (e.g. 
Chengalur, Salpeter & Terzian, 1993). This value is further 
justified by considering that, if M^ir > Mup in the above- 
specified luminosity range, then, the amount of dark matter 
locked in spiral galaxies fls = /^'/''fj ^vir 4>{.L) dL, with 
<j>{L) the galaxy luminosity function, would much exceed 
0.1, and would be unacceptable for the assumed cosmologi- 
cal model having fim = 0.3. 

We performed the fit to the data with the Xtot rnini- 
mization technique described in §3; the results are reported 
in Table 3 , Col.(l)...(5). In Figure 6 we compare the CDR 
and the NFW models: for 7 galaxies the NFW model is 
unacceptably worse than the CDR solution, whereas for 2 
objects (M-3-1042 and N7339) the goodness level of the two 
different fits is comparable, but the virial mass in the case 
of CDM is quite high: M,.„- ~ 2 x IO^^Mq. Moreover, the 
CDM models have extremely low value for the disk mass- 
to-light ratio (see Table 3 Col.5). The majority of the ob- 
jects have Md/Li < 0.05, in obvious disagreement with 



the spectro-photometric properties of spirals, that indicate 
mass-to-light ratios at least 10 times higher. 

The inadequacy of the CDM model for our sample 
galaxies is even more evident if one performs the fit after 
removing any constraint on virial mass. Indeed, (results in 
Table 3, from Col. 6 to 9), good fits are obtained only for 
very low values of the concentration parameter (cvir — 2) 
and for uncomfortably large virial velocities and masses 
(Kir ^ 600 - 800 km s'^Mvir ^ lO" - IO^^Mq). These 
results can be explained as effect of the attempt of the mini- 
mization routine to fit the V{r) oc r"'^ NFW velocity profile 
to data intrinsically linear in r. 



6 SUMMARY AND CONCLUSIONS 

We have performed the disk-halo decomposition for a well 
suited sample of 9 bulge-less disk galaxies with 100 < Vopt < 
170 km s"'^. These galaxies have a relevant amount of dark 
matter: the contribution of the luminous matter to the dy- 
namics is small and it can be properly taken into account. 
Moreover, the high spatial resolution of the available rota- 
tion curves allows us to obtain the separate dark and lumi- 
nous density profiles. We find that dark matter halos have a 
constant central density region whose size exceeds the stel- 
lar disk length-scale Rd- These halo profiles disagree with 
the cuspy density distributions typical of CDM halos (e.g. 
Navarro, Frenk & White, 1997; Kravtsov et al., 1998), which, 
therefore, fail to account for the actual DM velocity data. 
On the other hand, these halo velocities are well described in 
terms of the Burkert density profile, an empirical functional 
form whose two structure parameters (central density and 
core radius) are related through: po ~ fg^^^. 

This work is complementary to that of Salucci (2000) 
who derived for 140 objects of different luminosity V^, the 
logarithmic gradient of the halo velocity slope at Ropt- Re- 
minding that Vfi = and Vh — 1 mean an isothermal and a 
solid-body regime, respectively, let us stress that the result 
found: < 1 completely supports, at the optical edge, the 
results of the present paper obtained for a small sample of 
spirals over the entire stellar disk. 

Finally, the dark halo velocity linear rise from 0.25Rd to 
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^ 3Rd sets a serious upper limit to the dynamical relevance 
of CDM-like dark halos in spirals. Indeed, once we rule out 
a CDM halo, also the claim by Burkert and Silk (1997) of 
two dark halos, a MACHO dark halo with a CDR profile 
and a standard CDM halo, meets a difficulty. Indeed, in this 
case, in order to satisfy eq.(lO), the MACHO halo should 
account for > 95% of the dark mass inside Ropt- Then, it 
would dominate the dynamics out to well beyond 25 Rd- 
The CDM halo would then have dynamical importance in 
regions so external that its cosmological role itself would be 
in question. 
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